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We employ a kinetic-exchange tight-binding model to calculate the magnetic interaction and 
anisotropy energies of a pair of substitutional Mn atoms in GaAs as a function of their separation 
distance and direction. We find that the most energetically stable configuration is usually one 
in which the spins are ferromagnetically aligned along the vector connecting the Mn atoms. The 
ferromagnetic configuration is characterized by a splitting of the topmost unoccupied acceptor levels, 
which is visible in scanning tunneling microscope studies when the pair is close to the surface and 
is strongly dependent on pair orientation. The largest acceptor splittings occur when the Mn pair 
is oriented along the (110) symmetry direction, and the smallest when they are oriented along 
(100). We show explicitly that the acceptor splitting is not simply related to the effective exchange 
interaction between the Mn local moments. The exchange interaction constant is instead more 
directly related to the width of the distribution of all impurity levels - occupied and unoccupied. 
When the Mn pair is at the (110) GaAs surface, both acceptor splitting and effective exchange 
interaction are very small except for the smallest possible Mn separation. 



I. INTRODUCTION 

Experimental progress in the past 5 years has led 
to a large number of experimental studies of magnetic 
and nonmagnetic transition-metal impurities in semicon- 
ductors using advanced scanning tunneling microscope 
(STM) techniques.ElEI This effort was motivated in part 
by the hope that the high-resolution imaging and spa- 
tially resolved spectroscopic power of the STM could help 
in developing an accurate microscopic picture of dilute 
magnetic semiconductor (DMS) magnetism. In DMSs 
magnetic impurities provide local moments, which can 
couple to yield a collective ferromagnetic state. In the 
prototypical DMS, (Ga,Mn)As, the Mn impurities act 
as acceptors providing itinerant holes that can mediate 
long-range interactions between local moments. Among 
the open issues in DMS physics^ are the precise charac- 
ter of the hole states, the nature of the coupling of the 
holes with the local magnetic moments, and the prop- 
erties of the ensuing magnetic interaction between local 
moments. STM experiments performed recently^ ^ are 
playing a decisive role in clarifying some of these issues. 

In Ref. STM substitution techniques were used to 
incorporate individual Mn atoms into Ga sites in a GaAs 
(110) surface. Real-space spectroscopic measurements in 
the vicinity of an isolated Mn impurity revealed the pres- 
ence of a mid-gap resonance arising from a Mn induced 
acceptor state. High-resolution imaging showed that the 
acceptor wave function is strongly anisotropic with re- 
spect to the crystal axes of the host. When two Mn 
atoms were incorporated close to each other, two reso- 
nances appeared in the gap, split by approximately 0.5 
eV. The splitting was found to be strongly dependent 
on the Mn pair orientation with respect to the GaAs 
crystal structure and on Mn separation. A simple toy 
model, describing acceptor states coupled to the Mn ion 
local moments of the two impurities, suggested that a 



measurable splitting of the acceptor levels could only oc- 
cur if the two Mn local moments were ferromagnetically 
aligned.^ Since acceptor-level splitting is an observable 
indicator of ferromagnetic coupling, it seemed plausible 
that the dependence of this splitting on separation should 
be related at least qualitatively to the Mn-Mn exchange 
interaction. If so, the STM experiment could be used to 
measure exchange interactions between the Mn moments 
and test theories of this interaction. One of the purposes 
of the present study is to examine this relationship quan- 
titatively. 

Several of the experimental features uncovered in 
Ref. [4 could be qualitatively accounted for theoretically 
by a tight-binding model calculation for Mn in bulk GaAs 
presented in the same paper. A more thorough compari- 
son between experiments and theoretical modeling, both 
for Mn atoms in the bulk and Mn near the surface is 
nevertheless necessary to interpret the experiments, mo- 
tivating the present theoretical work. Here we consider 
a kinetic p-d exchange tight-binding model^^ in which 
the effective exchange interaction between the hole states 
and the local Mn moments arises from hybridization of 
the impurity d levels with p levels of the host. The model 
is solved numerically for large super-clusters containing 
up to 3200 atoms. This approach allows us to place the 
Mn pair either in bulk GaAs or on the (110) surface. 

Within this model, we study the electronic and mag- 
netic properties of Mn pairs in GaAs, assuming that the 
two local magnetic moments are collinear, having either 
parallel [ferromagnetic (FM)] or antiparallel [antiferro- 
magnetic (AFM)] relative orientation. One of our goals 
is to study the spin-orbit induced magnetic anisotropy 
energy of the system and see how this quantity is re- 
lated to the properties of the mid-gap acceptor states. 
In the FM configuration the magnetic moment tends to 
point along the direction of the pair, while in the AFM 
configuration there is typically a quasi-easy plane per- 
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pendicular to the pair direction. As in the isolated Mn 
impurity case previously studied}^ we find that the sum 
of the individual anisotropy energies of the top two unoc- 
cupied valence band orbit als, mirrors the total anisotropy 
of the system. This shows that the picture of Mn-Mn in- 
teractions mediated by valence band holes is valid, and 
simplifies the interpretation of our numerical results. We 
then consider the properties of the acceptor levels for the 
two possible relative orientations of the Mn magnetic mo- 
ments. In the FM configuration, which is generally found 
to be the ground state for most Mn pair orientations and 
Mn separations, the acceptor levels lie in a group just 
above the valence-band maximum and have a splitting 
that is enhanced by inter-ion hybridization. The group 
of six split levels can be viewed as a nascent version of 
the impurity band which forms in the bulk at small but 
finite Mn concentrations.'^ In particular, the two unoccu- 
pied acceptor levels (i.e. occupied by holes), which fully 
determine both magnetic interaction and anisotropy en- 
ergies have a finite splitting, which can be measured in 
experiment. For Mn pairs in bulk, we find that the ac- 
ceptor splitting varies strongly with pair orientation and 
Mn separation, in agreement with experiment and previ- 
ous calculations.^ The splitting is maximal for the most 
closely spaced Mn pair oriented along the (110) direction, 
where it is of the order of few hundred meV, and very 
small for pairs oriented along (100). For some Mn pairs, 
the wave functions of the two acceptor states have bond- 
ing and antibonding character, which is again consistent 
with experiment.^ ^ These results support the validity of 
the kinetic p-d exchange model. 

In our study the energy difference between FM and 
AFM moment orientations is related to partial occupa- 
tion of the acceptor impurity levels, which are split more 
widely in the FM configuration. When we calculate the 
exchange constant J for the Mn-Mn interaction, as the 
difference between the ground-state energies of the two 
configurations, we find that J is not in one-to-one corre- 
spondence with the splitting between the two unoccupied 
acceptor levels. In particular, the level splitting is much 
more anisotropic than J with respect to the Mn pair ori- 
entation. On the other hand, we find that estimates of 
the width of the distribution of all impurity levels - oc- 
cupied and unoccupied - are more directly related to J. 

Effective exchange interactions between Mn ions in 
bulk (Ga,Mn)As have been calculated by several groups 
using either ab initio methods or more phenomenolog- 
ical approaches. Comparison of our results with other 
estimates is not always straightforward, since we have 
only two Mn moments and the exchange interactions are 
not strictly pairwise. The order of magnitude of our J 
is nevertheless consistent with published results obtained 
from first-principles calculations, although values of J for 
specific Mn pair directions and separations may differ. 
However, it is well-known that DMSs are not accurately 
described by the local-density-approximation often used 
in first principles methods. The discrepancies could ei- 
ther be due to the shortcomings of our model or to the 



inaccuracies of ab initio calculations. 

We have also looked at how the Mn pairs interact when 
they are placed in a (110) GaAs surface. Typically, we 
find that the acceptor wave functions become highly lo- 
calized at the surface and produce states that are deep 
in the band gap. As a result, the long-range interactions 
are much weaker and antiferromagnetic alignment of Mn 
spins is more likely to occur. The magnetic anisotropy 
energies at the surface are an order of magnitude smaller 
than in bulk, and tend to produce quasi-easy planes at 
close distances, reverting to isolated Mn anisotropy land- 
scapes at larger separations. For Mn pairs in the (110) 
surface, the long-range behavior of the acceptor wave 
functions and acceptor splitting seems to agree less well 
with experiment than for Mn pairs in bulk GaAs. The 
tight-binding model that we use is of course less well 
justified when the Mn pair is in the surface. However, 
effects present in the experiment may contribute to this 
difference. In addition to an uncertainty in model pa- 
rameters at the surface due to band-bending shifts of the 
acceptors, a large overlap with continuum states due to 
Zn co-dopants in the sample can cause a more extended, 
bulk-like Mn acceptor wave function. Beyond the scope 
of the present paper, care should also be taken in simu- 
lating the change in the effective potential experienced by 
surface electrons due to the addition and removal process 
introduced by the STM. Our results clearly demonstrate 
that modeling STM studies of the surface of a system as 
complicated as (Ga,Mn)As is highly nontrivial and re- 
quires additional work. 

The paper is organized as follows. In Sec.|n|we review 
some theoretical aspects of the (Ga,Mn)As system and 
give a brief introduction to our tight-binding Hamilto- 
nian. We also elaborate on a toy model that gives an 
idea of the system behavior expected when the Mn mo- 
ments are parallel or antiparallel. The results for pairs 
of Mn along different symmetry directions in a fully pe- 
riodic bulk-like environment are presented in Sec. |III A[ 
and in Sec. IIIB for Mn pairs in the (110) GaAs surface. 
Finally, we summarize our conclusions in Sec. |IV[ 



II. THEORY 
A. Hamiltonian for Mn impurities in GaAs 

In this section we briefiy review the basic physics of 
Mn impurities in GaAs, which motivates the choice of the 
tight-binding Hamiltonian^^ used throughout the paper. 
For more details the reader is referred to Ref. [10 , where 
the same model was used to investigate the properties of 
individual Mn atoms in GaAs. 

In the neutral state, consisting of the Mn^+ ion on 
a Ga site with a weakly bound hole, the Mn has spin 
S = 5/2 and orbital moment L = OA^ This is in accor- 
dance with Hund's first rule that survives even when the 
Mn is embedded in the crystal. The Mn 3d^ electrons 
build up large, localized magnetic moments. Because 
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FIG. 1: Schematic drawing of the level dynamics leading to 
a polarized acceptor at the top of the valence band. The Mn 
3d^ up and down electrons are split by the atomic exchange 
in accordance with Hund's first rule. The levels are further 
split by the crystal field into a doublet of Cg and a triplet of 
t2p -symmetry. The -levels hybridize with nearest neighbor 
dangling bond As p-states, such that the levels at the top of 
the valence band with spin parallel to the Mn spin are pushed 
up in energy. Spin up is denoted by + and spin down by -. 



the Mn is missing the Ga 4p valence electron, it intro- 
duces a hole in the system, simultaneously acting as an 
acceptor and a source of magnetic moments. In bulk, 
the holes introduced by many Mn impurities form itin- 
erant carriers that mediate the ferromagnetic coupling 
between localized moments. Fig. [l] shows a schematic 
drawing of what happens when a Mn substitutes a Ga in 
the GaAs lattice. The Mn 3d^ up states (for example) 
that lodge in the valence band are exchange- split from 
the down states that end up in the conduction band. 
The crystal field imposes the tetrahedral host symmetry, 
resulting in a further split into a doublet of Cg symme- 
try and a triplet of t2g symmetry. The Cg states couple 
only weakly to the host. The t2g states on the other 
hand, hybridize with the As nearest neighbor dangling- 
bond p states at the top of the valence band, forming 
bonding and antibonding combinations. Effectively, the 
neighboring As p spins at the top of the valence band 
that are parallel to the Mn spin are shifted up in energy, 
whereas p spins that are antiparallel to the Mn spin are 
shifted down. Therefore, the main effect of a single sub- 
stitutional Mn with spin up, is to introduce three spin- 
polarized (" up" ) levels above the top of the valence band 
and three spin-polarized ("down") levels below. These 
six states are predominantly of As p character. The hole 
introduced by a single Mn will occupy the highest of the 
three antibonding states (indicated by the empty circle in 
Fig. [T]). In this way, the hybridization mechanism gives 
rise to an antiferromagnetic coupling between itinerant 
and localized spins by what is known as Zener's kinetic- 



exchange mechanism. ESHni account for this p-d in- 
direct exchange by introducing a classical Mn vector of 
magnitude 5/2, which couples to the neighboring As p 
states. There are no explicit d electrons in our model; 
they are accounted for only implicitly via the exchange 
interaction between the localized moments and the As p 
spins on neighboring sites. 

Our tight-binding kinetic exchange Hamiltonian takes 
the following form: 

ij,lili',a m n[m] 

i,IJ,/j,' ,aa' 

+ z^EE^H^+^co„. (1) 

ineoSr "Zf"^ fi-Rm\ 

The first term reproduces the band-structure of bulk 
GaA^IS] and contains the near- neighbor hopping an d on- 
site energies in terms of the Slater-Koster parameter Jl^'^ 
t^^^f for the s and p electrons of Ga and As. In (jlj) z and 
j are atomic indices, /i and fi' are orbital indices and a 
denotes spin. In simulating the (110) surface, buckling is 
accounted for by rescaling^^ the tight-binding param- 
eters and modifying the direction cosines appropriately. 

The second term implements the p-d exchange mecha- 
nism described above. It couples the unit spin vector of 
Mn atom m: (i^, to the p orbitals of the nearest neigh- 
bor As atoms n[m]: Sn = \ Y^^^a' «La^^^'^n7ra', where 
r is a vector of Pauli matrices. The value of Jpd has 
been inferred by Xheoiy'^^ and experiments^ to approx- 
imately Jpd = 1.5 eV. The three p-d hybridized levels, 
spin-polarized in the direction of the Mn moment (im 
are split from the three levels polarized in the opposite 
direction by an energy of order J^^- 

The third term in Eq. ([T]) accounts for SO interactions 
in an atomic aj^proximation with the renormalized spin- 
orbit splittings^^ Aca = 2AMn = 58 meV and Aas = 140 
meV. Spin-orbit interaction causes the band energy to 
depend not only on the relative angles between different 
Mn spin directions, but also on spin-orientations relative 
to the lattice. 

The fourth term of the Hamiltonian accounts for a 
long-range repulsive Coulomb interaction in the presence 
of a Mn impurity, which attracts a weakly bound hole 
and repels electrons. At the surface, we crudely account 
for the weaker dielectric screening by reducing the bulk 
dielectric constant Sj. in half. The Coulomb behavior at 
short distances is parametrized in Korr, which is a Mn 
central cell correction term used as a parameter to tune 
the Mn acceptor level to the experimental positionl^-^^^ 
It contains an on-site V^n = 1-0 eV (that is, acting on 
the Mn) and an off-site term VqE = 2.4 eV. The off- 
site Coulomb correction affects the nearest neighbor As 
atoms and together with Jpd controls the amount of p-d 
hybridization in the system and the range of the acceptor 
wave function. 
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FIG. 2: The six spin-polarized levels resulting from the p- 
d hybridization around a Mn impurity, shown in Fig. ^ The 
three levels with orbital index /i = —1, 0, +1 above the valence 
band edge, polarized in the direction of the Mn moment, are 
higher in energy than the three corresponding levels polarized 
in the opposite direction by an amount of the order of the 
effective exchange constant Jpd- Spin-orbit interactions (and 
surface effects when present) lift the orbital degeneracy of the 
like-spin levels and cause a small admixture of opposite spin 
character. The topmost up-spin level is the acceptor level 
introduced by the Mn. 
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FIG. 3: Toy model for two Mn atoms when the two local 
moments are aligned antiparallel - the AFM configuration. 
Each Mn site has six exchange-split itinerant levels that hy- 
bridize with the like-spin levels at the other Mn site. Note 
that like-spin levels have different energies at the two sites. 
The resulting hybridized levels, drawn in the middle, have 
their energy lowered or raised by ~ T^/JpcZ where T is the 
intersite hopping energy scale. This results in a net gain in 
energy because of the two holes introduced by the two Mn 



B. Toy model for Mn pairs in GaAs 

In this section we discuss a system of two Mn atoms 
in GaAs by means of a simple toy modeP that eluci- 
dates the basic mechanism responsible for the effective 
coupling between their localized spins. The properties 
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FIG. 4: Toy model for two Mn atoms when the two local mo- 
ments are aligned parallel - the FM configuration. Each Mn 
has six exchange- split itinerant levels that hybridize with the 
like-spin levels at the other Mn via hopping. Now, like-spin 
levels at the two sites are degenerate and give rise to bond- 
ing and antibonding states split by the hopping parameter 
2T which causes a "widening of the band". The resulting 
hybridized levels are drawn for the case when the hopping 
parameter T is smaller than the spin-orbit induced splitting 
^. Note that for the acceptor (hole) states above the valence 
band edge, we use the convention of calling "bonding" the 
state with the higher energy (subscript h — bonding and a — 
antibonding) . 
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FIG. 5: The same FM alignment of Mn spins as in Fig.|4] but 
for the case when the hopping parameter T is much larger 
than the spin-orbit induced splitting ^. The resulting impurity 
of bandwidth ^ 2T, is partially filled when two holes 
are present, resulting in a large energy gain with respect to 
the uncoupled T = system. 



of this simple model will be very useful in interpreting 
the results of our numerical calculations in Sec. IIIII 

For each Mn impurity in GaAs we will focus on the 
six spin-polarized levels at the top of the valence band, 
emerging from the p-d hybridization shown in Fig. [l] 
These six levels are shown explicitly in Fig. [2] for one 
particular orientation of the Mn moment. The three lev- 
els above the valence-band edge with orbital angular mo- 
mentum projection ja = 0,±1, are spin-polarized in the 
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direction of the Mn magnetic moment. Similarly, the cor- 
responding three levels below the valence band edge are 
polarized in the opposite direction. Spin-orbit coupling 
(and surface effects when the Mn is not in the bulk) lifts 
the orbital degeneracy of the like-spin states. We model 
this splitting by introducing the energies 



(2) 



where ^ represents the SO-induced splitting. Here (+) 
and (— ) stand for the levels above and below the valence- 
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band edge, respectively. Note that 
where Jpd is the effective exchange coupling used in 
Eq. ([T]). In the presence of spin-orbit interaction the spin 
is no longer a good quantum number, but we will assume 
that these states still have a predominant spin character, 
which is the same as when SO interaction is absent. 

When two Mn impurities are present and close to each 
other, the system can lower the total energy by allow- 
ing hopping between two single-particle states with the 
same spin, each centered around one of the two Mn. Two 
different situations arise depending on whether the rel- 
ative orientation of the two Mn moments is parallel or 
antiparallel. We assume that each state with 
a =t,4 at Mn site i will be coupled to the correspond- 
ing same spin state | j, /i, a) at the other Mn site j by 
an effective hopping parameter T^, which we take to be 
spin-independent. The single-particle Hamiltonian rep- 
resenting these two sets of itinerant spins coupled to the 
Mn local moment with hopping between sites is given by 
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where the upper signs in superscript refer to the FM and 
the lower signs to the AFM configuration.^^ The Hamil- 
tonian in E q. (p| corresponds to the Anderson-Hasegawa 
HamiltoniaiP^^ith the angle 6 between the spins equal 
to and tt, for the parallel and antiparallel configurations 
respectively. 

The Hamiltonian is immediately diagonalized by not- 
ing that the spin and orbital characters are good quan- 
tum numbers. The nature of the resulting hybridized 
levels will depend on the relative orientation of the two 
Mn moments: in the FM configuration the two sets of 
unperturbed equal-spin states (one at each Mn site) are 
degenerate and will form bonding and antibonding com- 
binations via hopping. By contrast, in the AFM con- 
figuration like-spins at different Mn sites have different 
energies and the hybridization will be reduced by a factor 

^ ^/i/ Jpd- 

In the AFM configuration the spectrum is doubly de- 
generate with energies 



>FM(±) _ AFM(±) 



jy^- (4) 



Fig. [3] shows a schematic view of the energy levels for the 
AFM configuration with a constant = T. In the pres- 
ence of holes, the total energy of the AFM state, ^^™, 
obtained by summing the energies of the occupied states, 
is lower than for the non- hybridized (T^ = 0) state. The 
maximum gain occurs for the "half-filled" system con- 
sisting of six holes. For the case of two Mn introduc- 
ing two holes (shown in Fig. |3|, the energy gain of the 
AFM configuration is approximately 2T'^ / J^^- This phe- 
nomenon is similar to the superexchange mechanism - 
arising within a one-band Hubbard model at half filling 
- which favors antiferromagnetic alignment of the itiner- 
ant spins on neighboring sites. Note that there are no 
d electrons present in our model and that the superex- 
change between the local moments is brought about by 
the kinetic exchange between the itinerant spins. When 
the two Mn atoms are close to each other, oppositely 
aligned Mn spins allow the wave functions to spread out, 
thus lowering their kinetic energy by hopping. 

In the FM configuration the spectrum is nondegenerate 



/M(+) _ (+) 
FM(-) 



(-) 



FM(+) 
FM(-) 



.(+) 
.(-) 



(5) 
(6) 



where £™ ^-^^ and £™ ^-^^ are bonding and antibonding 
states respectively. The acceptor (hole) states are more 
bound when they are further away from the valence band 
top, which is the start of the continuum for hole states. 
Out of the two acceptor states arising from the hybridiza- 
tion of the degenerate like-spin states of energy e^^\ we 
therefore assign^^ the label "bonding" to the one that 
occurs at higher energy, while the one with lower energy 
is denoted as "antibonding" [see Eq. [s] . 

The energy spectrum for the FM configuration is 
shown in Fig. [4] when the hopping parameter is smaller 
than the SO splitting, = T <C ^. The six accep- 
tor levels above the valence band edge form an "impu- 
rity band" with an associated bandwidth that increases 
with the splitting T. Note that for T <C ^, when two 
holes are present the total energy is the same as for 
the non- hybridized case, T = 0. An energy gain equal 
to — T occurs only when an odd number of holes are 
present. The opposite limit of a large hopping parame- 
ter, = T ^ ^, is shown in Fig.[5j Since only two thirds 
of the "impurity band" is filled there is a net energy gain 
^ — 2T due to the level splitting, which can stabilize the 
FM configuration against the AFM one. This mechanism 
that couples the itinerant spins ferromagnetically, corre- 
sponds to double exchange between the two sites. For 
a partially occupied impurity band the FM alignment 
tends to be energetically more favorable than the AFM 
when the widening of the band induced by hopping is 
large. The difference between the total energies of the 
two configurations defines an effective exchange constant 
J = (^^t"^^^ - Ej^)/2 representing a Heisenberg-like 
magnetic interaction, 1-L cx — J5'Mn(i) ' 5'Mn(i)7 between 
the local moments at the two Mn sites. 
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FIG. 6: Toy model for two Mn sites in the FM configuration 
when one hopping matrix element (T^^+i) is much larger than 
the other two and larger than the spin-orbit induced splitting 
^. Only the impurity levels above the valence-band edge are 
drawn. The half-occupied impurity band favors a FM align- 
ment over the AFM configuration and the non-hybridized 
(T/j^ = 0) system. The acceptor splitting is large and equal to 
Tu--\-i — T^=o+^. As we will see later (following the discussion 
10), this is the case that best describes the numerical 
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results of the microscopic Hamiltonian in Eq. ([T]). 



The actual hopping paths are of course more compli- 
cated than the ones shown in Figs. [SjSj In particular, 
for a hopping parameter depending strongly on the 
orbital character /i, we expect the resulting FM config- 
uration to be intermediate between the limiting cases of 
Fig. [4] and Fig. [5j In any case, the toy model predicts 
that the FM configuration will always be characterized 
by a splitting Aacc of the two acceptor states, related ei- 
ther to covalency between individual acceptor levels or 
to spin-orbit coupling. The splitting is noticeably absent 
in the AFM configuration, where the two acceptor states 
are quasi-degenerate, and is therefore a landmark of the 
FM state. As mentioned in the introduction, the STM 
experiments do measure a significant splitting between 
the two acceptor states, which is a strong indication that 
the Mn pair is coupled ferromagnetically. It is there- 
fore important to investigate if some kind of relationship 
exists between Aacc and J. This can be the case if, for 
example, the level structure of the FM configuration is of 
the type sketched in Fig. [6| Here a dominating hopping 
term Tn=+i gives rise to both a partially filled impurity 
bancP^ stabilizing the FM state, and a large acceptor 
splitting Aacc = ~ + ^, which will also be 

approximately related to the exchange energy gained by 
the FM configuration. As we will see later (following the 
discussion of Fig. 10), this is the case that best describes 



the numerical results of the microscopic Hamiltonian in 
Eq. 0. 

In the next Section we examine the properties of the 
six impurity levels when the Mn spins are parallel or an- 
tiparallel, and see how these states relate to the effective 
exchange constant J and the magnetic anisotropy energy. 
Note that the trace of the p-d exchange operator, summed 
over all valence band states is zero. If the impurity levels 
involve negligible conduction band character, it follows 



[010] [100] 

,(^ V^[110] 



[110] 
[001] 



\ I ^ 



<"fi>*V*.V(iio) 



FIG. 7: (Color online) Directions and separations of a Mn pair 
on the (110) surface. The two Mn atoms, replacing Ga atoms 
(bright green spheres), are marked (first Mn atom) and 1,2,3 
(second Mn atom) in correspondence with Fig. [s] below. Note 
that all equivalent symmetry directions (...) have been chosen 
in direct coordinates [...], such that all the Mn pairs appear 
in the (110) plane. 



that the sum of the energies of all valence band states 
is independent of Mn spin orientations.^^ Because of this 
property, exchange and interaction energies are expected 
to be accurately expressed in terms of the energies of the 
two empty valence band states. 

In the evaluation of J, we will also approximate the 
total energy for the FM and AFM configurations by sum- 
ming up the energies of the two unoccupied acceptor lev- 
els, or alternatively of the four topmost occupied levels. 
These results are compared the J obtained by the differ- 
ence between AFM and FM total energies. We will also 
consider other measures of the acceptor level structure 
broadening in Mn dimers, which might be more closely 
related to exchange interaction than the splitting of the 
top two levels. These measures include the splitting be- 
tween the mean of the four occupied and the mean of the 
two unoccupied levels, and the effective "bandwidth" of 
the six levels as obtained by calculating the standard de- 
viation from their mean value. 



III. RESULTS 
A. Mn-Mn interactions in bulk GaAs 

In this section we study numerically the properties of 
substitutional Mn pairs in bulk GaAs. As shown in Fig.[7| 
we consider Mn pairs oriented along different crystalline 
directions at various separations. The two Mn are em- 
bedded in a 3200 atom GaAs supercell with periodic 
boundary conditions in all directions, corresponding to 
a Mn fraction of 0.06%. We consider collinear magnetic 
configurations in which the Mn moments are either par- 
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allel or antiparallel. 



1. Magnetic anisotropy energy 

We begin by looking at the magnetic anisotropy of the 
system, which is defined as the dependence of the to- 
tal ground-state energy £^tot(^) on the direction of the 
Mn-pair magnetic moment, Ct. For the AFM configura- 
tion, Cl is the direction of the staggered moment. Graphi- 
cal representations of the magnetic anisotropy landscapes 
£^tot(^) on the unit sphere of the all possible directions 
for are presented in Fig. [s] for the FM configurations. 
Each panel (a)-(k) refers to a Mn pair with a particular 
separation and orientation in the crystal, according to 
the notation defined in Fig. [7[ln these figures, what we 
actually plot is 



£^anis(^) = ^tot(^) 



rrimin 
^tot 5 



(7) 



as a function of the magnetic moment direction Cl. Each 
point of the "anisotropy surface" is obtained by mea- 
suring the "distance" (£^tot(^) — ^tot") froi^i the center 
of the reference parallelepiped defined by the cubic axes 
[100], [010], and [001]. The center of the parallelepiped 
is also taken to be the center of the unit sphere of direc- 
tions. Here E^^^ is the minimum value of £^tot(^) upon 
varying Cl that occurs along one of the easy directions. 
The actual values of the high and the low energy bar- 
riers of these systems are shown in Fig. |9] In the FM 
configurations, the bistable easy directions are generally 
found to be parallel to Mn pair axis. The two excep- 
tions are the quasieasy planes perpendicular to the con- 
necting line, formed when the Mn are spaced at two and 
three lattice constants along (100) [see Fig. 8](f) and (g)]. 
For these two configurations, the magnetic hard direction 
has 12-15 meV higher energy. The qualitatively different 
landscapes signal that the interactions along the (100) 
direction differ from the other crystalline directions. For 
the shortest separation along (100), d = a = 5.65 A, a 
bistable easy axis parallel to the Mn-connecting line with 
a blocking barrier of around 11 meV is found. This means 
that a level crossing at the Fermi level occurs when the 
distance is increased from one to two lattice constants. 
Associated with this crossing is a change in the orbital 
character of the acceptors that results in a qualitatively 
different anisotropy. 

Focusing now on the (110) [with the Mn pair along 
[110], see Fig. [s] (b)-(d)], for which closely spaced As and 
Ga provide more direct hopping paths between the Mn, 
we see that there is a low barrier along the [001] axis 
and a high barrier along [110]. At the closest spacing, 
d = 0.7a, the anisotropy energy is very small with bar- 
riers of 1-3 meV [see Fig. [s] (b)]. This configuration is 
special, because the off-site Coulomb correction on the 
common nearest neighbor As between the Mn is additive, 
giving it a large on-site energy that reduces anisotropy. 
The (110)* pair in Fig. [s] (a) shows the effect of a non- 



(a)(110)S.0.71a 

[001], 
[100]- 
[110] 
[010] 



(b) (110)^^o.71a 




(e)jioo)^ (f)(ioo),.2a (g)(ioo),=3a 










(h) {2n)d=1.22a (i) {2n)d=2A5a 





(j) (211)^3.67a (k) (111)^^1 




73a 




FIG. 8: (Color online) Magnetic anisotropy in collinear vari- 
ation of parallel Mn spins on the unit sphere. Generally, an 
easy axis is found to be parallel to the connecting Mn line, 
with the exception of the (100)d=2a,3a- The (110)* shows the 
effect of non-additive Coulomb correction on the common As 
neighbor, yielding larger anisotropy barriers. In panel (a) the 
crystal axes have been marked and the same labels apply to 
all other panels. Also indicated is the (110) plane that con- 
tains all Mn pairs. The Mn are numbered in accordance with 
Fig.[7| The AFM configurations all have a higher total energy 
and the variations are inverted with respect to the FM ones. 
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1.0 1.5 2.0 2.5 3.0 3.5R/a 

FIG. 9: High and Low magnetic anisotropy barriers. In these 
graphs, zero corresponds to the minimum energy on the unit 
sphere. The anisotropy increases when the Coulomb correc- 
tion is taken to be non- additive for (110)* in the FM con- 
figuration (a). The (110)d=i.4i shows enhanced anisotropy 
barriers, with strong interactions along this symmetry direc- 
tion. (100) has a single barrier for the closest distance, but 
exhibits quasi-easy-planes (low barrier is zero) perpendicular 
to the connecting Mn line for longer distances. In the AFM 
configurations (b), the anisotropics gradually increase with 
distance as cancellations between the two acceptor levels de- 
crease. 



additive^^ Coulomb correction on the common As neigh- 
bor, where the low and high barriers have increased to 
9 and 14 meV. As we let the Mn move apart one step 
further along the (110) with distance d = 1.4a, the high 
barrier attains the maximum value of all considered con- 
figurations of 21 meV, and a low barrier of 12 meV. At 
the longest considered separation along (110), the low 
and high barriers are now 12 and 16 meV, respectively. 
The (211) Mn pairs follow a similar evolution of the high 
and low barrier, relative the connecting Mn line, now 
with a low barrier along the [110] direction that changes 
to the high barrier with increasing distance. Distances 
are longer than for the (110) pairs and therefore hy- 
bridization is weaker and the anisotropy energies drops 
by a few meV. At the larger separations along (211) the 
bistable minima are tilted away from the connecting line, 
an indication that the single Mn anisotropy is becom- 
ing comparable to the Mn-Mn interaction energy. Even 
with a supercell of 3200 atoms we cannot exclude pos- 
sible finite-size-effect contributions to the anisotropics. 
Nevertheless, we can get a good estimate of the barriers 
separating the generally bistable minima. The magnetic 
anisotropy in (Ga,Mn)As nanostructures, such as (Ga, 
Mn)As epilayers, is presently a topic of great interest. 



Understanding the microscopic mechanisms of mag- 
netic anisotropy in DMS is crucial in order be able to 
manipulate the magnetization vector by magnetic and 
electric fields. For a review of recent theoretical and ex- 
perimental works on magnetic anisotropy in (Ga, Mn)As 
see Ref. [31], where ferromagnetic samples with a rela- 
tively high Mn content of a few percent are examined. 
Comparison with our results, which instead pertain to 
the magnetic anisotropy of isolated Mn dimers, is not 
straightforward. Our results could on the other hand be 
directly compared with those of STM experiments using 
a spin-polarized magnetic tip^^ and an external magnetic 
field. 

As we show below, all the AFM pairs have a higher to- 
tal energy, which is expected on the basis of the heuristic 
toy model considerations. The AFM anisotropy land- 
scapes are qualitatively inverted with respect to the FM 
ones, i.e., the bistable minima along the Mn dimer axis 
in the FM variation are replaced by quasieasy planes per- 
pendicular to the dimer axis in the AFM variation. This 
inversion of hard and easy directions between FM and 
AFM can be understood by considering the FM configu- 
ration in the easy direction, where the spin-orbit splitting 
between the highest occupied and the lowest unoccupied 
levels is large, resulting in a lower total energy. In the 
corresponding AFM configuration, the two levels are of 
opposing spin character, such that the splitting is smaller 
and the total energy instead higher. In Fig. [9](b), we see 
that the closer the Mn are, the smaller the anisotropy en- 
ergies tend to be. This figure gives us an idea of the range 
of strongly interacting Mn. Weaker interactions at larger 
distances in the AFM configuration tends to increase the 
anisotropy, as the effectively spin-polarized region in the 
lattice around each Mn increases. 



2. Character of acceptors for Mn pairs in bulk GaAs 

The anisotropy of the embedded Mn dimers is accu- 
rately reflected in the anisotropy of the two acceptor lev- 
els and their associated wave functions. In this section 
we therefore perform a detailed analysis of various prop- 
erties of the acceptor levels, such as their splitting and 
LDOS - quantities that are directly accessible by STM 
spectroscopy. 

From now on we use the convention to denote with 
E-s ,£-2 ,£-1 ,£o the energies of four highest occupied 
states, while ei and 62 denote the energies of the two 
acceptor levels. When mixing with the conduction band 
is negligible, the anisotropy energy landscape can be ex- 
tracted from the the two acceptor levels using 

2 

E^^U^) = -J2{sm - ef% (8) 

i=l 

where £^anis(^) is defined in Eq. ^ and the upper limit 
in the sum can be generalized to a greater number of ac- 
ceptor levels. Equation ([8| follows from the fact that the 
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trace of the p-d exchange operator in Eq. ([T]) summed 
over ah valence band states including the acceptors is 
zero. In the FM variation of the Mn spins, the high- 
energy acceptor level, £2{^)^ varies very little and in a 
qualitatively opposite manner to the low-energy accep- 
tor, ei{Cl). Its effect is therefore to reduce the much 
higher anisotropy coming from the lower acceptor level. 
In the AFM configuration small acceptor level variations 
partially cancel to a low total anisotropy below 1.5 lattice 
constants, and add up to a larger total anisotropy above 
this value. 

SO interaction mixes spin components, such that the 
eigenstates are no longer of definite spin character. As a 
result, the acceptor states above the valence band edge, 
which without SO coupling have the same spin charac- 
ter as the corresponding localized Mn moment ("spin- 
up"), acquire a small component in the opposite direc- 
tion ("spin-down"). In the FM cases, this results in a 
small spin down character of the acceptor levels, that in- 
creases with Mn distance primarily for the low acceptor, 
as it moves closer to the valence band. In the AFM varia- 
tions, the two acceptor levels are now very close in energy 
and their spin character can vary by a large amount on 
the unit sphere due to spin-orbit interaction. 



In Fig. 10 we plot the energies of the four highest oc- 
cupied and the two unoccupied acceptor levels for all 
studied configurations. We consider spherical averages 
of these six energies over all possible moment directions, 
denoted by Si , i = —3 , . . . ,2. 

In the FM configuration [Fig. 10 (a)] we find that there 



is a clear splitting between the three lowest and highest 
levels, in agreement with the results of the toy model 
shown in Fig. [5] The fact that three highest levels are 
also spread over a sizable energy range indicates that, at 
least for the shortest Mn separation, we are in the regime 
described by Fig. [6] of the toy model. Note in particular 
the large splitting of the two acceptor levels, as expected 
for aligned Mn spins. The level splitting and the overall 
width of the impurity band^^ decreases with increasing 
Mn separation. Note also the strong dependence of the 
acceptor splitting on the orientation of the Mn pair. In 
particular, the FM acceptor splittings tend to be largest 
along the (110) direction. It is also clear that the (100) 
Mn pair behaves differently. Even if the splitting is large 
between the three lower and three upper levels, the split- 
ting of the two acceptor levels is consistently smaller than 
for the other orientations. 

In the AFM configuration the structure of the six en- 
ergy levels is quite different and its salient features are 
nicely captured by the toy model result of Fig. [Sj The 
levels are always essentially doubly degenerate; the split- 
ting between doublets can be large for the shortest Mn 
separation and but it decreases quickly with distance and 
becomes of the order of the expected SO-induced split- 
ting. In contrast to the FM case, there is no visible split- 
ting between the acceptor levels. The only exception is 
the special case of the (110) for the shortest Mn separa- 
tion, where the six impurity levels abruptly drop down 
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FIG. 10: The four highest occupied (filled circles) and the 
two lowest unoccupied (empty circles) eigenstates as obtained 
by taking the spherical average. In the FM configuration 
(a) the splitting between the three upper and lower impurity 
levels decreases with distance, but the split between the two 
acceptor levels follows a more complex pattern, depending on 
pair orientation. In the AFM configuration (b) levels bunch 
up but are still split from each other by weaker hybridization. 
The (110)* has been shifted slightly to the right in order not 
to obscure the graph. 



towards the valence continuum as the off-site Coulomb 
correction on the common As neighbor is decreased. Fi- 
nally note that the (100) pair sticks out with a more 
dense set of levels also in the AFM configuration 

We proceed to discuss the Local Density Of States 
(LDOS) of the two acceptor levels - a property that can 
be probed by STM spectroscopy-l^l^] plots of the LDOS 
of the two acceptor levels in the (110) plane containing 
the two Mn, are shown in Fig. [TT] for the FM and in 



Fig.[T2]for the AFM configurations. The LDOS plots are 
generated^^ from the tight-binding model calculations as 
in Refs. [10 33 . In these plots the Mn spins are pointing 
in the magnetic easy directions, indicated by the arrow 
in the top left corner of each panel. The left column in 
each figure shows the LDOS of the lower energy acceptor 
state, which is closer to the valence band maximum; the 
right column refers to the acceptor with higher energy. 
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FIG. 11: (Color online) LDOS of the lower-energy (left col- 
umn) and the higher-energy (right column) acceptor level in 
the easy direction (indicated by white arrow in top left cor- 
ners) for the FM configurations. Each row corresponds to a 
given crystal orientation and separation for of the Mn pair. 



FIG. 12: (Color onhne) LDOS of the lower-energy (left col- 
umn) and the higher-energy (right column) acceptor level in 
the easy direction for the AFM configurations. 
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We first consider tiie (110) pairs in the FM configu- 
ration (see Fig. 11). For the special case where the Mn 
atoms are separated by a common nearest neighbor As, 
the results change significantly depending on whether the 
off-site Coulomb correction is additive {{llO)d=o.7ia) or 
non- additive ((110)5=o.7ia)- ^^^^ effect of a 

non-additive Coulomb correction on the common As is to 
delocalize the acceptor wave functions. A more delocal- 
ized wave function generates larger anisotropy energies, 
as shown in Fig. [9j At the shortest Mn separation, both 
acceptor wave functions have bonding character^ with the 
maximum spectral weight on the common As neighbor 
located between them. This situation is captured by 
the toy model, where the large hybridization occurring 
at small Mn separation gives rise to the level structure 
shown in Figs. |5] and [6} with both acceptors being of the 
bonding type. As the Mn ions move apart along (110), 
the lower acceptor state develops more antibonding char- 
acter with a significant decrease of the spectral weight 
between the two Mn sites, whereas the upper acceptor re- 
mains in a bonding-like state. Within our toy model, this 
implies that the two acceptor states correspond to bond- 
ing and antibonding states arising from the hybridization 
of degenerate levels in energy, spin and orbital character, 
as shown in Fig. [4j 

In agreement with our results, the STM experiments 
for the (110) pair find clear evidence that the two Mn- 
induced acceptor states have bonding and antibonding 
character, with the bonding state occurring at higher 
energies.^ ^ An antibonding character for the lower en- 
ergy state is observed experimentally also for the Mn pair 
with the shortest separation.^ While this does not seem 
to be case in our bulk (110)ci=o.7ia calculations (both 
states being essentially "bonding"), we do find that the 
value of the maximum LDOS on the As in between the 
Mn's is around twice as large for the upper acceptor wave 
function. 

For the (100) pairs the upper acceptor exhibits a more 
localized signature, whereas the lower acceptor is more 
extended because it is closer to the valence continuum. 
Bonding and antibonding formations cannot be clearly 
seen, indicating that the hybridization between the two 
holes is much weaker along the (100). This conclusion is 
also supported by the very small acceptor-energy split- 
ting, as shown in Fig. [2l] In the remaining directions 
(211) and (111) (rows 8-11), bonding and antibonding 
characters are visible, with more and less spectral weight 
in the region between the Mn, respectively. 



> 



In the AFM LDOS (see Fig. 12) there are no bonding 
or antibonding patterns for the two acceptor states, in 
agreement with the results of the toy model. Rather, a 
spatially symmetric separation seems to occur for some 
of the Mn pairs, in particular for larger distances along 
the (100), (211) and (111) directions. For these pairs the 
upper and lower acceptor acquire opposite and rather 
definite spin characters in the easy direction. In this 
case spin up is located on one site and spin down on the 
other, whereas for (110) the two acceptors are of mixed 
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FIG. 13: Acceptor level splittings as a function of pair orien- 
tation and distance, (a) shows a comparative plot of current 
and theoretical and experimental results from Ref. [4]. In (b) 
and (c) the FM and AFM results for all considered pairs are 
shown. 



spin character. 



S. Acceptor splitting vs. effective exchange constant J 

In this Section we take a closer look at the acceptor 
energy splitting Aacc(^) = ^2 — ^i, focusing in particular 
on its spherical average Aacc = S2 — Si. We will also 
compute the effective exchange constant J{tt) and see if 
a relationship can be found between these two important 
quantities. 



Fig. 



13 



J shows the average acceptor splitting Aacc in the 
FM and AFM configurations compared to the experimen- 
tal and theoretical values reported in Ref. [4^. Our cal- 
culated splittings match the previous calculation^^, and 
follow the trend of the measured splittings at the (110) 
surface. Fig. 13 (b) reveals that the splitting is largest 
for Mn pairs along (110) directions, and very small for 
(100) Mn pairs. The (211) and (111) Mn pairs exhibit 
smaller splittings than (110), which has more direct hop- 
ping paths. Notice that there is some correlation between 
the splittings and the high barriers in Fig. |9](a); higher 
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FIG. 14: Effective exchange couplings. A ferromagnetic cou- 
pling is always favored as all spherical averages J > 0. 



anisotropy barriers correspond to a larger splitting. 
Fig. 13 (a) shows that our tight-binding model as well 
based on a similar approach, sys- 




as a previous theorj/^ 
tematically underestimates the acceptor splitting mea- 
sured experimentally. This discrepancy could be due to 
inherent limited accuracy of the tight-binding method. 
However, it is also quite possible that part of the splitting 
seen experimentally is of Coulombic origin arising when 
electrons tunnel into the two different acceptor levels. In 
Fig. [13] (c) we see that the AFM splittings are very small 
compared to the FM ones, and are typically of the order 
of a few meV. As already noted above, the only exception 
is the (110)* pair. Reducing the Coulomb parameter on 
the common As neighbor dramatically lowers the accep- 
tor energies [see Fig. (b )], produces more extended 
wave functions (see Fig. |12|), and increases mixing with 
valence band states. 

We can now consider the effective exchange energy de- 
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FIG. 15: Effective exchange couplings for a system with two 
Mn and one hole. Adding one electron to the system, gives an 
effective exchange that is related to the acceptor splittings. 



is totally absent in the plot for J. A similar dip for the 
{100)d=2a pair is hinted but much less pronounced in the 
J plot. 

At this point it is useful to consider other quanti- 
ties that can shed light on the relationship between 
the effective exchange interaction and the acceptor lev- 
els. We first consider the effective exchange constant 
for a Mn pair when only one hole is present or, equiv- 
alently, when the lower acceptor is occupied by an elec- 
tron. The toy model results in Figs. [3] and |4] can help 
intuition and justify the rational of this choice: when 
an extra electron is added to these electronic configu- 
rations, only the topmost acceptor level is unoccupied 
in the FM state. In this case the acceptor splitting 
should be directly related to the energy gain stabilizing 
a FM state over the AFM state. The effect of occupy- 



fined as J{n) = [E^^(^ (Q) - E[^^ {n)]/[2- where ing the lowest acceptor is shown in Fig. [Tsl where we 



the factors of 5/2 from the Mn spin magnitudes is ab- 
sorbed in J. A positive (negative) J implies that the 
Mn-Mn interactions are FM (AFM). As a result of the 
spin-orbit interaction, J{Cl) is an anisotropic quantity, 
that is, it depends on Cl. 

The effective J as obtained by taking the spherical av- 
erages is shown in Fig. [l4j where it can be seen that the 
FM configuration is always the most stable one, J > 0. 
Only for separations greater than two lattice constants, 
do the SO-induced fluctuations in J become comparable 
to the average value. Even then, ferromagnetism is al- 
ways favored. The maximal value of J generally occurs 
when the Mn spins are pointing along the pair axis. 

The comparison between the acceptor splitting A™ 
for the FM state [Fig. [is] (b)] and the exchange constant 
J [Fig. 14 (a)] shows that both quantities decay rapidly 
with Mn separation. However, there are noticeable differ- 
ences between them, out of which two are most obvious. 
First of all, the relative value of J for the (110)c^=o.7ia 
pair compared to all the other values is much larger than 
the corresponding value of A™ . Secondly, A™ displays 
a much less monotonic decrease with Mn separation than 
J. In particular, the large dip for the {100)d^a,2a pairs 



plot Jihole = m 



tot,lhole 



tot,lhole 



]/[2-(5/2p], in which 



^AFM,FM 
^tot,lhole 



is the spherical average of the total ground 
state energy for a system with one Mn pair and only 
one hole, namely with one extra electron added. We can 
indeed see that now the dependence of Jihoie on the crys- 
tal orientation and spacing of the Mn pair is qualitatively 
much more similar to the A™ of Fig. 13 (b), displaying 
the same large dips for the (100) pairs. 



In Fig. 16 we plot the spherical averages of the effective 
J as estimated by the FM/AFM difference of the sums 
of the 4 occupied impurity levels and the 2 unoccupied 
levels. Comparing with the "exact" J calculated by total 
energy difference of the full system of electrons (Fig.[l4]), 
we see that both approximations underestimate the value 
of J for the most closely spaced (110) Mn pair. For the 
rest of the points, using the four highest occupied levels 
overestimates the value of J approximately by a factor 
of two. Surprisingly, taking the difference of the sums of 
the two acceptors in the FM and the AFM configurations, 
gives a very good estimate of the "exact" J for the rest of 
the points in the plot of Fig, 14 (a). In order to compute 



the effective J for the two acceptors, one must invert the 
sign, because 2 unoccupied levels are used. 



13 



40 



30 



> 
S 



x20 

O 



10 



o 
o 

o i 



Oo 



O using 2 unoccupied 
• using 4 occupied 



O 



rsi 



o 



oo o 



o 



o 



0.5 1.0 1.5 2.0 2.5 3.0 3.5 R/a 

FIG. 16: The effective exchange constant as estimated by the 
FM/AFM difference of the 4 occupied and the 2 unoccupied 
acceptor levels. 
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and Fig. 17 (b) shows the splitting between the mean of 
the four occupied impurity levels and the mean of the 
two unoccupied levels, 



6 = 



i=l 



(10) 



In the FM configuration these can both be seen as a mea- 
sure of the effective exchange interaction strength. A 



comparison of Figs. 17 and 14 reveals that they both do 
quite well in reproducing the correct trend in J, with the 
exception of the most closely spaced (110) pair. These 
observations support the validity of the simplified argu- 
ments used in the toy model. 

In conclusion, our analysis of the acceptor splitting 
versus effective exchange constant shows that these two 
quantities are certainly related but not in a direct quan- 
titative way. In particular, the acceptor splitting is much 
more anisotropic than J as a function of the Mn pair ori- 
entation. Other quantities such as w and 6 [see Eqs. (|9| 
and (10)] involving all six impurity levels seems have a 
much better correspondence with J. It is interesting to 
ask whether or not the four occupied impurity levels are 
also accessible by STM. One might naively expect that 
when the applied bias is reversed, electrons tunneling out 
of the occupied impurity levels should appear as sharp 
features in the differential conductance spectrum. The 
experiments of Ref. [4 do not show any clear sign of these 
levels. However, a more recent stud}/^^ on individual Mn 
impurities in GaAs indicates that a weak feature due to 
tunneling out of one the three t2g levels can indeed be 
seen below the top of the valence band. Presumably this 
is due to a hole- hole interaction effect: the final state has 
two electrons missing near the Mn impurity. This feature 
is very hard to see, as it is masked by electrons tunneling 
out of the valence band continuum. 



4' Comparison with ah initio estimates of the exchange 
coupling J 



FIG. 17: Six impurity level bandwidth (a), and splitting be- 
tween the mean of the four occupied and the mean of the two 
unoccupied impurity levels (b). 



In our discussion of the toy model, we have seen that 
the emergence of a stable FM configuration for the Mn 
pair is brought about by the "widening" of a partially 
occupied cluster of impurity levels caused by hopping. It 
is therefore instructive to compare some measures of the 
"bandwidth" of the six impurity levels with J. Fig. 17 
(a) shows the FM effective bandwidth of the six impurity 
levels. 



w = 



E 



IT. 



(9) 



In this section we compare our results for the effec- 
tive exchange J with those of previous studies. Since the 
pair exchange interaction between impurities is a cru- 
cial quantity in the field of DMS, there have been many 
theoretical studies of this quantity, mainly based on first- 
principles calculations. Here we concentrate our atten- 
tion on a few issues that have emerged from our theoret- 
ical approach: (i) the relatively large value for the (110) 
pair with shortest Mn spacing; (ii) the non-monotonic de- 
cay of J with Mn separation and its strongly anisotropic 
character with respect to pair orientation; (iii) the com- 
parison with acceptor splitting, with particular reference 
to the discrepancies between J and Aacc found in our 
model, e.g., for the (100) pairs. A few caveats are nec- 
essary before comparing our results with other estimates 
of J that have appeared in the literature. Firstly, most 
of the published estimates are for much higher Mn con- 
centrations, which strongly affects the value of J and its 
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anisotropic properties. Secondly, all first-principle cal- 
culations suffer from the well-known limitations of DFT 
when applied to semiconductors. In particular, the es- 
timates of J depend on the version of DFT used, e.g. 
LSDA vs. GGA. When the GGA+U approach is used 
to account for electronic correlation effects, the value of 
the parameter [/, which has as strong influence on the 
results, can only be determined indirectly by comparing 
with experiment. 

In Ref. [35 the electronic structure of (Ga,Mn)As 
was calculated from first principles including disorder 
via the coherent-potential approximation. The magnetic 
force theorem and one-electron Greens functions then al- 
low mapping onto an effective Heisenberg Hamiltonian. 
Their results for the exchange interaction strength re- 
veal a strong dependence on the doping concentration x. 
In particular, the effective J of the nearest (110) pair 
is highly sensitive to doping and increases dramatically 
with decreasing Mn concentration, ranging from J = 19 
meV for x 0.08 to J = 55 meV for x 0.001. The 
lowest concentration of 0.1% Mn, agrees very well with 
our estimate. Note that our calculations correspond to 
a slightly smaller effective Mn concentration x = 0.0006. 
In Ref. |9 the authors employed a self-consistent LSDA 
atomic-spheres approximation^^ and their result of the 
effective J = 55 meV for the {H^)d=Q.7ia pair agrees 
remarkably well with our values, despite the small super- 
cell corresponding io x = 0.009. With spin-orbit interac- 
tion this value drops to 48 meV. In Ref. the authors 
employ the GGA+U method within the Projector Aug- 
mented Wave (PAW) ah initio approach, treating U as 
a parameter. For (110)c^=o.7i and x = 0.03 they observe 
a decreasing J as a function of [/, with J = 34 meV for 
U = and J = 22 meV for U = 6 eV. As the Mn concen- 
tration is reduced from x = 0.03 to x = 0.008 at /7 = 0, 
an increase in J of approximately 10 meV is observed. A 
larger U parameter causes the d levels to emerge deeper 
in the valence band, the acceptor wave function becomes 
more delocalized and p-d exchange decreased. Photoe- 
mission experiments^^ indicate that the t2^-states should 
be approximately 4 eV below the valence band maximum, 
which corresponds to U = 3 — 4 eV. Ref. [37] shows that 
the calculated J depends on the extension of the hole 
wave function as per chosen U. 

We now turn to the (100)c^=a pair, which has a very 
low experimental acceptor splitting (see Fig. 13). In con- 



trast to our res ults (see F ig. 14), ab initio calculations 
generally predictP^'^^I^^'^ a dip in the curve of J vs Mn 
separation occurring for this pair, with a lower value than 
the two following points, (211)^^=1. 2ia and {110)d=iAia- 
Similar but smaller dips are also found for the other (100) 
pairs.— However, as mentioned above, all these calcula- 
tions show a significant dependence on impurity concen- 
tration. For example, Ref. [35] finds that the value of 
J for the (100) d^a pair steadily increases from approx- 
imately J = 1 meV for X = 0.1, to J = 7 meV for 
X = 0.02. So do the next two points for (211) and (110) 
which end up at approximately 10 meV and 11 meV for 
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FIG. 18: The effect on J when increasing the Mn concen- 
tration from X = 0.0006 (shown in Fig. [m]) to x = 0.0017 
(filled circles). The empty circles show how J changes when 
increases x further for selected points. 



X = 0.02. When = and x = 0.008, the GGA -1- 
U approach finds J = 6 meV, well below the values for 
the next two pairs (211)^=i.2ia and {llO)d=iAia, J = 14 
meV and J = 9 meV respectively, which are also com- 
parable to our results. Upon increasing U to the physi- 
cally reasonable value of 6 eV, for both (211)^^=1. 2ia and 
(110)c^=i.4ia J = 7 meV, while J = 3.5 meV for the 
(100)d=a pair. 

For even larger values of U the difference in J between 
the (100)c^=a and (211)^^=1. 2ia pairs decreases further; 
however, the (100) (^=a is still consistently lower. The 
GGA+U values of J for {110)d^iAia and {'211)d^i.2ia^ do 
not agree with the experimental acceptor splitting trend 
for sound values of U. A similar trend for GGA calcu- 
lations is reported in Refs. [4^41. The LDA result 
for 1% Mn including disorder in the coherent-potential 
approximation agree reasonably well with our results for 
(100)rf=„ {2n)d=i.2ia and (110)^=i.4ia with J = 8.2, 9.1 
and 7.9 meV. The (100)c^=a is still consistently lower for 
larger distances and the (110)c^=o.7ia is relatively low at 
J = 18 meV, indicating a sensitivity of J on the chosen 
ab initio method. 

In summary, with the caveats mentioned above, the 
general trend of the ab initio J at reduced concentrations 
shows some qualitative agreement with STM experimen- 
tal acceptor splittings at shortest distances d = 0.71a, 
d = a. For larger Mn separations the agreement is not 
so good. As corresponding theoretical estimates of the 
splittings are not reported in the literature, it is hard 
to conclude whether or not a relationship between the 
acceptor splittings and J is present in the ab initio cal- 
culations. Concerning the comparison between the ab 
initio values of J and our results, we can see that al- 
though the overall magnitudes agree reasonably well, the 
qualitative trend differs since at short Mn separations we 
find a monotonic decrease of J with distance. 

Given the strong dependence of J on Mn concentra- 
tion and the fact that the reported ab initio results are 
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obtained for larger concentrations, it is interesting to in- 
vestigate how our value of J changes when we increase 



(a) (llO)^^o.na 



(b) (110>^^1.41a 



the Mn concentration. In Fig. 18 the effect of decreas- 
ing the supercell size such that the Mn concentration 
changes from x = 0.0006 to x = 0.0017 is shown. The 
(100)(i=a J is now closer to the (211)(i=i.2ia value. In- 
creasing the concentration further causes the two pairs 
to obtain equal J, but the (100)(i=a value never drops 
below the {HV)d^i.2ia value, indicating that there is a 
fundamental difference with respect to ah initio. 

Generally, ah initio predicts a qualitative trend that 
agrees better with our effective J obtained by occupying 
the lower acceptor (see Fig. 15). The (100)c^=a highest 
occupied level and lower acceptor are quasi-degenerate, 
with a gap that varies between 14-46 meV. In ah initio 
methods, a common technique used to speed up evalu- 
ation of k-point sums, is to introduce a fractional occu- 
pation of the unoccupied levels controlled by an occupa- 
tion smearing parameter. At the end of the calculation, 
the limit of zero smearing is taken to calculate the to- 
tal ground state energy. It is still possible that the end 
result depends on the choice of smearing parameter. If 
this is the case, the effect can cause the J for the quasi- 
degenerate (100) (i^a to decrease. 



B. Mn-Mn interactions in the (110) GaAs surface 



So far we have studied Mn pairs embedded in bulk 
GaAs. In this last Section we consider the experimentally 
more relevant situation where the two Mn atom substi- 
tute two Ga atoms on the (110) surface. We consider 
again a 3200-atom supercell, but this time we apply pe- 
riodic only in the two directions in the plane of the (110) 
surface. This corresponds to a 38 x 38 A^ surface and 
a supercell cluster that has 20 atomic layers along the 
surface normal separating the two (110) surfaces, such 
that surface-surface interactions are negligible. The loss 
of coordination and hybridization with the surface states 
will cause the impurity levels to appear very deep in the 
gap. The depth is subject to experimental uncertainty 
due to a band-bending^ effect when imaging a semi- 
conductor surface with a metal tip. In our model, the 
depth is very sensitive to the off-site Coulomb correc- 
tion, V^ff, and we use this parameter to reproduce the 
acceptor level of a single Mn on the (110) surface at the 
experimentally observed position^ at 850 meV. The pa- 
rameter VoE has therefore been reduced from 2.4 eV to 
1.57 eV. 

We proceed again to study the magnetic anisotropy 
energy of the systems, the properties of the mid-gap ac- 
ceptor states, the effective exchange interaction and its 
connection with the acceptor splitting. 




(e) {mT=0.7,a 



(f) (100),"^ 



AFM 




FIG. 19: (Color online) The magnetic anisotropy landscapes 
for Mn pairs in the (110) surface. The transparent plane 
indicates the (110) surface, (a)-(d) show the distinct types 
of landscapes in the FM configurations. Only along the 
(110)d=o. 71a, 1.41a in (a) and (b), and (100)d=a in (c), do the 
landscapes differ from the result of weakly interacting pairs 
in (d), which is just the sum of two isolated Mn anisotropics. 
Similarly, in the AFM configurations, only (110)d=o.7ia (e) 
and [100) d^a (f) deviate qualitatively from the landscape (d). 



1. Magnetic anisotropy energy 



Fig. 19 shows the qualitatively distinct anisotropy 



landscapes for Mn pairs on the surface and Fig. 20 the 
magnitude of the barriers. The surface geometry com- 
pletely dominates the anisotropy and for most separa- 
tions the anisotropy landscape is qualitatively the same 
as for a single Mn in the surfac^^^ [exemplified in Fig. 19 
(d)]. This type of anisotropy landscape occurring for sep- 
arations larger than ^1.5 lattice constants is an indica- 
tion that the acceptor state hybridization is so weak that 
the resulting anisotropics are essentially the sum of two 
independent Mn atom anisotropics. Only for the closest 
pairs do different anisotropics appear. In the FM con- 
figurations only (110)d^o.7ia, (110)d=i.4ia and (100)d=a 
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FIG. 20: The high and low barriers for the FM (a) and the 
AFM (b) configurations for pairs at the (110) surface. The 
anisotropy for the closest pairs exhibit a reduced anisotropy 
energy for separations greater than two lattice constants, all 
pairs have a characteristic barrier of around 2 meV, indicative 
of very weakly interacting Mn. 



[Fig. [19] (a)- (c)] show qualitatively different landscapes - 
tilted quasi-easy planes with the hard direction approx- 
imately along [111]. Similarly for the AFM configura- 
tions, only (110)d=o.7ia and (100)d=a [Fig. [l9] (e)-(f)] dif- 
fer from the single Mn type landscape in Fig. IT9| (d). 
In (100)^=^ the easy axis and the low barrier are inter- 
changed between the AFM and FM, but the hard direc- 
tion remains along the [111]. 

Fig. [20] shows that anisotropy energies are very small, 
typically one order of magnitude smaller than for the 
fully periodic systems. In both FM and AFM configu- 
rations the interactions between Mn for distances below 
two lattice constants, tend to reduce the anisotropy heav- 
ily. At larger distances, where all Mn pairs produce the 
[111] easy axis type landscape in Fig. [19] (d), the single 
barriers are around 2 meV. This limit of weakly interact- 
ing Mn is also reflected in the difference in barriers (see 
20) between the FM and AFM, which becomes negligi- 
ble for weakly hybridized pairs above 1.5a. As we will 
show, total energy differences are also small, causing a 
very small effective exchange. 



Character of acceptor levels for Mn pairs in a (110) 
surface 



Fig. 21 shows the group of six impurity levels for the 
Mn pairs in the surface. We see that for both FM and 
AFM configurations, two of the occupied levels are split 
around 0.5 eV below the others. For the closest (110) 
pair the acceptors are so high in energy that the upper 
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FIG. 21: The four highest occupied (filled circles) and the two 
acceptor levels (empty circles) in the FM (a) and AFM (b) 
configuration for Mn pairs in the (110) surface. Loss of coor- 
dination and hybridization with surface states cause impurity 
levels deep in the gap. 



acceptor has in fact crossed the first conduction band 
state. This is associated with a very large splitting of 
the acceptors. In (110)^ the off-site Coulomb has been 
reduced even further (to 0.5 eV) in order to also make the 
upper acceptor appear in the gap. The overall spectrum 
for the AFM and FM configuration looks similar, indi- 
cating that interactions are very weak and the itinerant 
spin wave functions are very localized at the surface. In 
the STM experimenlP Zn dopants can give rise to reso- 
nant tunneling between conduction states and the accep- 
tor level, thought to be responsible for the negative dip 
seen in the curve of the tunneling conductance vs. bias 
voltage-f^The coupling to conduction-like states can lead 
to a more extended, bulk-like acceptor wave function seen 
at the (110) surface. As a first approximation our treat- 
ment of the surface should reveal some of the relevant 
properties. Our results indicate that the interactions in 
the surface are much weaker than in bulk, and magnetic 
anisotropics are one order of magnitude smaller. This is 
related to acceptor wave functions that are much more 
localized at the surface. 

So far there are relatively few papers that attempt to 
simulate STM images of Mn in the (110) surface .^^^ The 
acceptor LDOS resulting from our tight-binding model 



for a few representative pairs are shown in Fig. 22 for 
the FM and Fig. [23] for the AFM configurations. The 
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FIG. 22: (Color online) LDOS for FM aligned Mn pairs on 
the (110) surface. The inset in (110)^ shows the level above 
the low acceptor when the upper acceptor has crossed into 
the conduction band. The images are for the Mn spin in the 
easy direction, indicated in the top left corner of each panel 
(note that [111] is coming out of the (110) surface). 
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FIG. 23: (Color onhne) Examples of LDOS for AFM aligned 
pairs Mn in the (110) surface. 



spectral weight in the core regions of the Mn are much 
higher than in bulk. As an example, the upper accep- 
tor for the (110) 0^71^ has around 50% spectral weight on 
the As in-between the Mn, to be compared with 17% in 
bulk, which is still high relative other bulk configurations. 
Wave functions for other surface pairs have a much higher 
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FIG. 24: The splittings of the acceptor levels for Mn pairs in 
the (110) surface. 



maximum spectral weight than in bulk, with values in 
the typical range 20%-30%. Overall, the spectral weight 
maximum for a given pair in bulk is approximately an 
order of magnitude smaller than the corresponding pair 
at the surface, with bulk wave functions that are much 
more spread out in the lattice. 

As shown in the first two rows of Fig. |22j some 
bonding/antibonding characteristics in the acceptor 
wave functions can be seen for the closest pairs 
(110)(i^o. 71a, 1.41a the FM configuration, which exhibit 
magnetic quasi-easy planes. Both acceptor wave func- 
tions for the {110)d=2.i2a Mn pair exhibit a hint of an- 
tibonding character. Note that this pair is characterized 
by the magnetic anisotropy landscape of a single Mn im- 
purity. 

For {100) d=a each acceptor level is occupying both Mn 
sites, but a separation occurs for {100)d=2a,d=3a such that 
each acceptor has more pronounced LDOS on one of the 
sites. The same association of acceptors spatially bound 
to one site, appears for the other directions - a pattern 
indicative of weak hybridization. Finally the (21 1)^^=1. 22a 
has some hint of bonding/antibonding pattern, but not 
as pronounced as for bulk. 



In the LDOS for the AFM configurations (see Fig. 23), 
the (110) pairs show no spatial separation between the 
acceptors and no bonding/antibonding behavior. The 
(211)c/=i.22a has upper and lower acceptor wave functions 
that are localized to the separate Mn sites, whereas the 
rest of the pairs are similar to the bulk counterparts (see 
Fig. 12), but with a more localized LDOS signature. 



3. Acceptor splitting vs. effective exchange constant J 

Finally, we have calculated the effective exchange con- 
stant for the Mn pairs in the surface, which we compare 
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FIG. 25: The effective exchange constant for Mn pairs in the 
(110) surface. The effective J falls off much faster at the 
surface, due to the highly localized acceptor wave functions. 
The (211)(i=i.22a is showing antiferromagnetic behavior. 



with the calculated acceptor splittings. The FM split- 
tings of the acceptor levels is shown in Fig. 24 (a). The 



very large splitting of 600 meV for the closest Mn pair 
seen in experiment (see Fig. 13) is now comparable to the 
calculated one of around 500 meV, when both acceptors 
are in the gap. The calculated splitting is small for the 
nearest (100). It then increases to around 100 meV for 
the (211), but after that splittings are very small and less 
than 50 meV, with the exception of the farthest (110). 
In the AFM configurations, splittings are much larger at 
the surface compared to bulk. In particular, the closest 
(211) has a large splitting of around 70 meV. This is ex- 
pected, since a finite splitting in the AFM configuration 
is possible when inversion symmetry is broken, an effect 
that is clearly enhanced by the presence of the surface. 

The FM surface splittings do not match the long range 
experimental trend as well as the bulk calculations. In 
our calculations, we have accounted for the buckling of 
the GaAs (110) surface, but not the detailed relaxation of 
Mn environment. It is also not certain that the effective 
p-d exchange model holds so well for the surface. Not 
only the long range and short range Coulomb may differ 
at the surface, but the exchange splitting is also subject 
to change. This is because the p-d hybridization scheme 
in Fig. [l] will be altered by the surface states. 

For the nearest 
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The surface J are shown in Fig. 
(110) pair, J is strongly ferromagnetic and has a value 
similar to the corresponding pair in bulk. The (100) pair 
has a very low J of around half an meV and could easily 
become antiferromagnetic by a change in Coulomb pa- 
rameters. The SO-induced anisotropy of J on the unit 
sphere of magnetic moment directions is very weak, just 
fractions of an meV, due to the low anisotropy energies 
at the surface. Interestingly the (211) pair, which has 
a relatively large acceptor splitting of several tens meV 
in the AFM configuration (see Fig. [24|), is antiferromag- 
netic. As mentioned above, this is an effect entirely due 
to the surface, where a significant acceptor splitting can 
exist even for antiferromagnetically aligned Mn spins and 



generate an energy gain that can compete with the FM 
configuration. The rest of the Mn pairs shown in Fig. 25 
have a very small J, indicating that the Mn local mo- 
ments interact very weakly due to the high degree of lo- 
calization of the acceptor wave functions at the surface. 
The oscillatory behavior is indicative of Friedel oscilla- 
tions with a much shorter period than in bulk, which 
is connected to the strong localization of the acceptor 
wave functions at the surface. Finally, it is interesting 
to note that the magnitude of J (disregarding the sign), 
follows the trend of the experimentally observed split- 
tings, despite the fact that the calculated ones do not 
clearly do so. Disregarding the single point (110)c^^i.4ia 
in Fig. 24 (a), the experimental trend but with smaller 
magnitudes is reproduced. It is also interesting to note 
the large AFM splitting for (211)^=1. 22a in Fig. [24] (a), 
which makes it difficult to distinguish from the FM con- 
figuration. There are however noticeable differences be- 
tween the (211)^=1. 22a FM and AFM LDOS (see Figs. 
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and 23). The FM acceptors each have spectral weight on 



both Mn sites, but the AFM acceptors are each associ- 
ated with one site. Given the large splitting between the 
AFM acceptors, this observation could perhaps be tested 
experimentally. Overall, it would appear that the exper- 
imental results for the (110) surface are intermediary to 
our surface and bulk calculations. 



IV. CONCLUSIONS AND OUTLOOK 

In this paper we have investigated substitutional Mn 
pairs in GaAs, using a microscopic tight-binding model 
solved for large GaAs clusters. The model accounts 
realistically for spin-orbit interactions and includes the 
kinetic-exchange p-d coupling between the Mn local mag- 
netic moment and the spin of valence band electrons. 
One goal of this work was to assess whether or not this 
model can reproduce the main features of coupled accep- 
tor states observed in recent STM experiments,^ ^ and 
in this way further elucidate the nature of the effective 
exchange interaction between the two Mn magnetic mo- 
ments. When the two Mn atoms are located in the bulk, 
we find that the ground state of the system is generally 
one in which the two Mn magnetic moments are coupled 
ferromagnetically and aligned along the direction con- 
necting the Mn atoms. In the ferromagnetic configura- 
tion, the two topmost acceptor states are split, and their 
energy separation depends strongly on the Mn pair ori- 
entation with respect to the GaAs crystal structure and 
decreases with Mn separation. In particular, the largest 
splitting {r^ 300 meV) is found for the Mn pair oriented 
along the (110) direction. For separations larger than 
three lattice constants, the splitting is below 20 meV. 
The splitting is the result of hybridization of degener- 
ate acceptor states with identical spin character, giving 
rise to bonding and antibonding states. Our calculated 
acceptor splittings are typically a factor of two smaller 
than the experimental values. Coulomb correlation ef- 
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fects not included in our tight-binding approach could 
partly be responsible for the larger splittings. Analysis 
of the LDOS for the two acceptor states shows that these 
are mainly concentrated in the surroundings of the Mn 
atoms; for intermediate separations, when the splitting is 
sizable, the higher state and lower state are of the bond- 
ing and anti-bonding type respectively. These results 
are all consistent with experiment,^ ^ indicating that the 
TB model correctly describes the electronic properties of 
the coupled acceptor states in GaAs associated with the 
magnetic impurities. 

An important question that we have addressed here 
is whether or not the acceptor splitting, which is acces- 
sible by STM measurements, is directly related to the 
effective exchange coupling between the Mn local mo- 
ments. The present work demonstrates that, at least 
within the model considered here, this relationship is not 
very sharp: although both quantities decrease with sep- 
aration and are typically anisotropic, the acceptor split- 
ting displays a non-monotonic behavior for a few pair 
orientations which is not found in the exchange coupling. 
In our model, the dependence of the latter on Mn sepa- 
ration and pair orientation is better represented by the 
overall bandwidth^^ of the six impurity levels present in 
the gap (see Fig. [l7). In particular. Our calculations 
demonstrate that exchange interactions between two Mn 
ions in (Ga,Mn)As are closely related to the splitting be- 
tween the two acceptor levels that are occupied by holes 
and the four lower energy acceptor levels that are occu- 
pied by electrons. This splitting is not readily measured 
by STM experiments which are strongly influenced by 
Coulomb interaction energies when the number of holes 
bound to the acceptor complex is increased above two by 
removing electrons from the system. It is possible that 
infrared spectroscopy of transitions between iso-charge 
levels of the acceptor complex could be more effective. 

These results support the idea that the stability and 
strength of the ferromagnetic coupling involves the hy- 
bridization of all impurity levels and not only the two 
topmost ones occupied by holes. Our estimate of the 
exchange coupling is in qualitative agreement with re- 
sults obtained from ab initio calculations, with a few no- 
ticeable discrepancies analyzed in detailed in Sec. |III A"4l 
which are most likely due to the limitations of both meth- 
ods and deserve further investigation. 

When the two Mn atoms are inserted in the (110) sur- 
face as in experiment, we typically find much deeper 



acceptor states with strongly localized wave functions 
around the impurities. Weaker hybridization causes the 
acceptor splitting to decay more rapidly with separation, 
although in a manner qualitatively similar to bulk. With 
the exception of the most closely spaced Mn pair, the 
surface acceptor splittings seem to agree less well with 
experiment than our bulk results. The exchange con- 
stant is likewise very small except at the shortest sepa- 
rations and can even change sign, favoring an antiferro- 
magnetic coupling. These results are a clear indication 
that our semiphenomenological microscopic model, de- 
rived mainly from bulk (Mn,Ga) As properties, is not able 
to quantitatively capture the complexity of the surface 
states. This is an area that could benefit from additional 
work. 

The model investigated here makes predictions for the 
spin-orbit induced magnetic anisotropy for pairs of Mn 
atoms in GaAs. Some of these predictions could be tested 
in STM experiments by applying an external field. When 
the Mn pair is in a bulk environment, there is a uniax- 
ial magnetic anisotropy along the axis of the pair with 
anisotropy barriers of 10-20 meV. When the pair is in a 
(110) surface, the behavior of the anisotropy energy land- 
scape is for most pairs similar in character and magnitude 
to that of a single Mn atom in the (110) surface. 

In the present approach the local moment degrees of 
freedom have been treated as classical variables. The 
study of the quantum spin dynamics of two interacting 
Mn in GaAs is an important and interesting subject both 
theoretically and in view of future experiments. Work in 
this direction, based on Berry-phase quantization of the 
semi-classical local moments considered here, is under- 
way. 
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